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We study the dynamic structure factor of two coherent bright solitary waves in attractive Bose- 
Einstein condensates confined in a harmonic trap. We demonstrate that the wave function of the 
two solitary waves with a fixed relative phase shows interference in momentum space. The fringes 
are shifted depending on the values of the phase. This momentum interference can be extracted from 
the dynamic structure factor of the system using stimulated two-photon Bragg scattering. Thus our 
method provides a way to measure the relative phase directly. 



I. INTRODUCTION 

' Bose-Einstein condensates (BECs) with attractive two- 
body interaction are dynamically unstable in a trap if 
the number of condensed atoms reaches a critical value 

■ [11, 0, i, H ■ Generally, hefty BECs will collapse following 
a switch of the s-wave scattering length from a positive 
value (or zero) to a negative one via the use of a Feshbach 
resonance Q . However it was found that in such case 
the number of remnant condensate atoms can be above 
the critical number [1, 0, S This was explained by 
the generation of brig ht solitary waves with repulsive in- 
teractions [Tol [ill Il2| . These repulsive interactions have 
been used to indirectly infer a repulsive relative phase 
n/2 < 9 < 37r/2 [l^, [HI, [i2| between neighbor solitons 
[3, [lj| . Origin of such phases could be quantum mechan- 
ical phase fluctuations of the bosonic field operators [l3| 
due to the nonlinear modulation instability. 

On the other hand, quantum and thermal excitations, 
which cannot be taken into account by simple mean- 
field theory, play a crucial role in the collapse of BECs 
[1, [3) d) [^ • These effects have been studied by effective 

I quantum field method [ij, [3, [H, [13] ■ An interesting re- 
sult arising from these studies is the absence of 7r-phase 
relation as revealed recently by [l^. Briefly, a relative 
phase of tt is introduced in mean-field theory between 
neighbor bright solitons [|. Usually one expects a tt 
relative phase to support repulsive interaction between 
solitons. However, there are indications that the repul- 

I sive interaction can come from quantum (thermal) noise 
rather than the repulsive phase [1^. In addition one 
also doesnot get a fixed 7r-phase relation within mean- 
field theory if BECs is developed from a nonuniform ini- 
tial state fly]. It is therefore interesting to ask what a 
meaningful phase would entail. A solid starting point 
would be mean field theory, in which 7r-phase is assumed 
[lOl [Til . [T3 ]. CoUisional dynamics of bright solitons in 
attractive BECs confined in a trap is usually based on it 
[T^. Yet other phase relation is argued in [l3,[l^ well. 
Consequently ascertaining the phase becomes important. 

In this paper, we propose a method to measure the rel- 
ative phase directly rather than to infer it indirectly from 
soliton interactions. In the experiment [rj, two bright 



solitary waves are separated in coordinate space. This 
makes it hard to measure the phase by overlapping the 
density [20| . However interference appears if the wave 
function is transformed into momentum space [2l[. This 
indicates that we can extract the relative phase from the 
dynamic structure factor. The dynamic structure factor 
provides important information about the spectrum of 
collective excitations and momentum distribution [22l ]. 
Stimulated two-photon Bragg scattering has been used 
to measure the dynamic structure factor at high momen- 
tum transfer and provided high resolution and sensitivity 
[23| . Measuring phase in momentum space is proposed 
by Pitaevskii et al. [2l|, who suggest to measure the 
coherence of two spatially separated BECs confined in 
a double-well potential in this way. This method can 
be applied for the problems we consider here as our two 
bright solitons are well separated in the trap [rj. We 
calculate the dynamic structure factors for two solitons 
with TT-phase and in-phase {9 = 0) using the impulse ap- 
proximation [24,, 25\ . Our results show that the dynamic 
structure factors have distinguishable peaks correspond- 
ing to these phase relations. 



The paper is organized as follows. In sec. II, we show 
that the momentum density exhibits interference due to 
the relative phase between two bright solitary waves in a 
trapped BEC. The phase carried by solitary waves is suc- 
cessfully recovered by a maximal response of the dynamic 
structure factor. Only when 9 = appears a single peak, 
which represents a large occupation of scattered atoms 
with the same energy as the photon. In contrast, two 
peaks would be the sign of 6* = tt. In sec. Ill, the dy- 
namic structure factor is discussed for a case when two 
solitary waves move symmetrically in a trap with a small 
velocity. This is the general case as they are trapped in 
a harmonic potential and interact repulsively. Although 
the shapes of the dynamic structure factor are modified 
by the nonzero velocity, they provide distinguishable in- 
formation of the concrete values of the phase also in this 
case. In sec. IV, we conclude the paper. 
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II. DYNAMIC STRUCTURE FACTOR FOR 
SOLITARY WAVES IN AN ATTRACTIVE EEC 

The Gross-Pitaevskii equation (GPE) is very effective 
in providing both quahtative and quantitative predic- 
tions for static and dynamical properties of an attrac- 
tive BEG O, H HI- We consider N weakly interact- 
ing bosons trapped in a quasi-one-dimensional potential 
Vtrap — m{ujlx'^ + ujpp'^)/2, where p'^ — y'^ + and 
LUp ^ LUx are the transverse and the longitudinal trap 
frequencies accordingly. The dynamics of the condensate 
is governed by the GPE 
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V' +Vt 



trap 
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here Uq — ANnti^aQ/m is the effective coupling constant 
representing two-body interaction strength, ao and m are 
s-wave scattering length and mass of atoms. The wave 
function of the condensate is normalized as J \ ip\'^dr — 1. 

The GPE is widely accepted for studying dynamics of 
interacting BECs at temperature T « OK. Explicit ana- 
lytical solutions are only available in a few situations. For 
instance, bright soliton solutions can be derived in free 
space using the inverse-scattering method 28|. Vortex 
solutions in rotatin g B ECs are also available for analyt- 
ical treatment pol. l30l|. It is hard to obtain the exact 
soliton solution in the presence of harmonic trap. In this 
paper, we consider the soliton solution of a BEG with a 
Gaussian form as it is confined in a harmonic trap. The 
wave function of a BEG confined in an elongated axially 
symmetric harmonic trap is given by [10, 31] 



^o(r) = 
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where li is the width of condensate in the z-direction and 
tti = y^h/moji, i = x,y,z are the harmonic lengths cor- 
respondingly. Note that the Gaussian wave function is 
close to the exact soliton solution when Lj; is comparable 
with Ux [13, HH. With this Gaussian ansatz, the energy 
of each atom is obtained readily from the GPE 



N ^ 



En 



Eint 



C/n 



(3) 



here we have used Eq and Eint denoting the zero-point 
energy in a harmonic trap and two-body interaction en- 
ergy. 

A homogeneous condensate with attractive interaction 
in free space is unstable while the zero-point kinetic en- 
ergy can prevent the collapse of a trapped condensate 
with the finite number of particles [2? 



at given trap frequencies and s-wave scattering length, 
a condensate containing a finite number of particles is 
stabilized by the balance of kinetic energy and negative 
interaction energy. On the other hand, the interaction 
energy is limited by a small value too. This will turn out 
to be helpful in calculating the dynamic structure factor. 

In the presence of two bright solitary waves, the wave 
function is taken as V(r) = ipo{r — ri) -I- e'^Volr — ^^2) 
with relative phase 9. Because of repulsive interactions 
between two solitary waves, they will be spatially sepa- 
rated in a harmonic trap [y, 0|. This makes it difficult 
to measure the phase in coordinate space, where overlap- 
ping density profiles could induce a fringe pattern [20l |. 
This obstacle can be overcome if we transform the wave 
function into momentum space, where two spatially sep- 
arated condensates will exhibit interference |21j. And 
the momentum distribution can be measured by stim- 
ulated two-photon Bragg scattering [23'| giving the dy- 
namic structure factor of the system. 

At large momentum transfer q, we can use the impulse 
approximation to calculate the dynamic structure factor 



S{q,E) = J 



= dpS[E- 
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where q and E are momentum and energy transfered 
from the laser pulse to the condensate and n{p) is the 
momentum density of the BEG 



n{p) = ((/.(p)V(p)) 



(6) 



where (pip) is momentum distribution obtained by 
Fourier transformation of wave function V-'(r) 



0(P) 



(7) 



In order to obtain the dynamic structure factor, at 
first we need to know the momentum density distribu- 
tion. We use the Gaussian ^ to calculate the mo- 
mentum distribution. For two solitary waves moving 
in the x-direction, the wave function of the system is 
-0(r) = i)o{x + 1/2, p) + e^^ipo{x - 1/2, p). I is the dis- 
tance between them. The corresponding wave function 
in momentum space is 
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where 0o(p) reads 
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where c is a constant of order 1 and a 
the mean harmonic length. It implies that on one hand. 



For the above wave function, the momentum density 
exhibits interference 



n(p) = 2[1 + cos(6' - pxl/h)]nQ{p) 



(10) 
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FIG. 1: (Color online.) Dynamic structure factor shows 
fringes for different phase. Solid, 6 — 0, dashed, 8 — n; 
dotted, a single condensate. Other parameters are, I = 
18^171, Ix = 7iim,q — 23ft/im~^. Ifere v — E/2-Kh and 
Vr = Er/2-nh. 
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FIG. 2; (Color online.) Same with Fig. HI but with vr x 10% 
relative phase deviation from 9 = Q,-n. Dynamic structure 
factors shift with respect to Er- Thick curves show the result 
of 6* + TT X 10% and the thin ones show 9 — n x 10%. 



where no{p) ~ |0(p)p. The period of the interference 
pattern is Dp — 2'Kh/l. Obviously momentum den- 
sity distribution depends on the phase 0. For example, 
TT-ipx = 0,Py,Pz) = at 9 = Pxl/fi + TT but has a maxi- 
mum at 6 = pxl/h. This implies that the relative phase 
9 displaces the momentum interference. 

In the following, we study the dynamic structure factor 
of the system by applying Bragg spectroscopy [23| . Here 
interaction effects on the relative shift of the peak are 
neghgible (Eq. ([3])). Weak interaction is particular useful 
because it makes using of im pulse ap proximation feasible 
at high momentum transfer |22l . |24| . 

By applying a Bragg pulse in the x-direction, we obtain 
the dynamic structure factor from Eq. ([5|) 

fn f 

Siq,E) = — J dpydp^n{px,Py,Pz) (H) 
with px = m{E — /2m)/q. Substituting Eq. pUj) into 



the above equation, we get 

S{q, E) = 2[1 + cos(0 - pxllh)]Ss{q, E) (12) 
Here Ss{q, E) of a single solitary wave is 

mix ^i'^^'^ it,'^ 



Ssiq,E) 
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(13) 



Eq. (jl3p has a Gaussian form. Like the momentum den- 
sity, Eq. (fT2|) also shows an interference pattern. The 
period of cosine function is 
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At a given distance / and optical momentum q, the 
relative phase 9 will shift the interference pattern of the 
dynamic structure factor. Consequently, we can measure 
the relative phase from the fringe of the dynamic struc- 
ture factor. In the experiment [7], the maximal distance 
of two solitons is about ISfim and the width of each soli- 
ton is about 7fim. Since the size Ix of an attractive con- 
densate is small, the wave length of the laser field must 
be far smaller than Ix ■ Here we assume the momentum of 
the laser field is q — 23fi./ir7i~^. In Fig. ^ we show the 
dynamic structure factor for two different phase relations 
9 = 0,n. When 9 = 0, S{q, E) is close to that of a single 
condensate (22l . [25| . In contrast to the maximum located 
at E = Er {Er = q^/2m) for 9 = 0, S{q,E) for 6* = tt 
it is zero there due to the interference. But two peaks 
emerge symmetrically around energy E = Er- These 
results show that the dynamic structure factor exhibits 
contrasting interferences for = and 9 = tt. 

As we introduced previously, the relative phase plays 
a crucial role in collapsed condensates 0, [lO, [13] ■ In a 
Bose gas with attractive two-body interaction, only a fi- 
nite number of atoms can be condensed Q. This is true 
for a single condensate. However, in a fragmented con- 
densate, for example, several bright solitary waves inter- 
acting repulsively support a bigger system. In this case, 
the repulsive interaction among neighboring peaks comes 
from the special relative phase 7r/2 < 9 < 3tt/3 accord- 
ing to mean field theory [T^, [ll|, . On the other hand, 
when including quantum noise, truncated Wigner quan- 
tum field simulations give 9 = rather than 9 = tt [l8l |. 
Since the experiments did not directly monitor the phase, 
a measurement proposal based on the dynamic structure 
factor is established to distinguish the relative phase. 

We also study the dynamic structure factor if 9 is 
slightly deviating from 9 = and tt. The deviation is 
10% X TT around them. For two bright solitons, a relative 
phase in the entire range 7r/2 < 9 < 3tt/2 can stabilize 
the system. Such a phase deviation could be possible 
in the experiment. Consequently, we'd like to see how 
a slight phase deviation affects momentum interference 
and whether or not we could distinguish the phase. The 
results are reported in Fig. ([2|). The dynamic structure 
factor remains similar except a slight displacement of the 
peak position. For phase shifts around 9 = 0, it only 
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FIG. 3; (Color online.) Dynamic structure factor for a pair 
of moving solitary waves. Solid, 9 = 0, dashed, 9 = n. Other 
parameters are, I = 13fj,m,Vx ~ 0.1mm/s,lx = 9,^ra,q = 



moves the maximum along the energy axis. The dynamic 
structme factors are not symmetric about E—Er — any 
more for phase shift about 6 — n. A maximum appears 
in the lower energy part for a shift 9 — (1 — 10%)x7r 
while other smaller peaks appear in the high energy part. 
Vice versa for 9 = {1 + 10%) x tt. Despite these shifts, 
the resulting momentum interference could give clearly 
distinguishable fringes for in-phase and 7r-phase. These 
studies demonstrate the robustness of the dynamic struc- 
ture factor against phase shift around 9 = t: and 9 = 0. 



III. DYNAMIC STRUCTURE FACTOR FOR 
MOVING SOLITARY WAVES 

Thus far we considered the situation when the two soli- 
tary waves have zero velocity and are well separated. 
However, the observed solitary waves move in the trap 
and the velocity cannot be always zero. Actually, the 
above results assume that the two solitary waves are 
separated by the maximal or minimal distance. In ex- 
periment two solitons move symmetrically about the 
trap center. The frequencies of the oscillation depend on 
the external harmonic trap frequencies. In this section 
we discuss the case when the two solitary waves have a 
small velocity oscillating in trap. The wave function for 
a moving solitary wave is 



-0fc(r) = il^oii 



ikx/h 



(15) 



where Vc is the center of mass. This is a moving wave 
packet with momentum — fc. Transforming it into mo- 
mentum space, we get 



c(P) 
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FIG. 4: (Color online.) Same with Fig. (O but with tt x 10% 
relative phase deviation from 9 = 0,7r. Dynamic structure 
factors shift with respect to Er. Thick curves show the result 
oi 9 + TV X 10% and the thin ones show 9 — tt x 10%. Other 
parameters are the same as in Fig. Q. 



where I = 2xc- Hereafter, we assume there are two waves 
oscillating symmetrically in the x-direction with momen- 
tum k = ±mv with respect to the trap center. The dy- 
namic structure factor for the two moving solitary waves 
can be obtained as 



Skiq,E) = —FiUkipx) 

q 

where nk{px) is the momentum density of Eq. JT 
the a;-directions and 



(17) 



E 



1 -I- exp 



-2 cos( 



Ipx , ( 2kpJl 
_)exp(^^^ 



(18) 



characterizes the interference of the dynamic structure 
factor. 

For a rough estimate, the maximal velocity for a par- 
ticle in a harmonic trap is Vmax ~ ujx^x, where Ax is 
the distance to the trap center. In the experiment 
Ax ~ 9^771, LUx = 2tt X 6.8Hz. One can get the maxi- 
mal velocity if it can reach the bottom of the potential, 
Vmax ~ 0.39mms~^. The velocity should be less than 
Vmax because of the separation of two solitary waves with 
repulsive interactions. Considering this point, the dy- 
namic structure factor with nonzero velocity is plotted 
in Fig. ([3]). We can see that the interferences are differ- 
ent from the zero velocity cases (Fig. ([T])). The maximal 
momentum distribution at E = Er is lowered for 9 = 
and the two peaks are bigger for 9 = tt relatively. But 
one can still distinguish the two situations with phase 
or TT. 

Again, we have calculated the situation when the phase 
deviates from 9 = 0{tt). These results are shown in Fig. 
(|3|), a slight shift of the maximal response will not obvi- 
ously reduce the contrast of the dynamic structure factor. 
Thus one can get the phase relation from measured data. 
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IV. CONCLUSION 

In this paper, we have discussed the dynamic struc- 
ture factor for two bright sohtary waves with different 
relative phases in attractive BECs. In mean field theory, 
the phase 6 = ir protects the condensate from collapse by 
providing a repulsive force between solitary waves. How- 
ever, both in-phase and 7r-phase relations can achieve this 
once quantum effects are considered [13, [HI [IE Ell ■ 
Our study shows that momentum interference is modi- 
fied by the phase. The dynamic structure factor displays 
distinguishable fringes at 6* = and tt respectively. At 
the proper photon energy and momentum, it has a sin- 
gle peak located at E — Er for 9 = while there are 
double peaks for 9 = tt. The distinguishable fringes are 



robust even when the solitary waves have a small veloc- 
ity. 10% phase deviations from 6* = and 9 = n are also 
discussed. Our result offers possible way to directly mea- 
sure the relative phase between two neig hboring solitary 
waves in attractive BECs [1, 0, [3, H 113]. An exper- 
iment using Bragg scattering could provide more detail 
information on the relative phase. 
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